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Topological superconducting states in monolayer FeSe/SrTiOa 
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The monolayer FeSe with a thickness of one unit cell grown on a single-crystal SrTiOa substrate (FeSe/STO) 
exhibits striking high-temperature superconductivity with transition temperature Tc over 65K reported by recent 
experimental measurements. In this work, through analyzing the distinctive electronic structure, and providing 
systematic classification of the pairing symmetry , we find that both s-and p-wave pairing with odd parity give 
rise to topological superconducting states in monolayer FeSe, and the exotic properties of s-wave topological 
superconducting states have close relations with the unique non-symmorphic lattice structure which induces the 
orbital-momentum locking. Our results indicate that the monolayer FeSe could be in the topological nontrivial 
s-wave superconducting states if the relevant effective pairing interactions are dominant in comparison with 
other candidates. 

PACS numbers: 74.70.Xa,74.78-w,74.20.Rp 


I. INTRODUCTION 


Topological superconductors [1-4] and iron-based 
superconductors [5] have been research focuses of condensed 
matter physics in recent years. Topological superconductors 
have a full pairing gap in the bulk and gapless surface or 
edge Andreev bound states known as Majorana fermions. Re¬ 
cent scanning tunneling microscopy/spectroscopy (STM/S) 
measurements observed a robust zero-energy bound state at 
randomly distributed interstitial excess Fe sites in supercon¬ 
ducting Fe(Te,Se), and the behavior of zero-energy bound 
state resembles the Majorana fermion[ 6 ]. Theoretically, one 
possible scenario accounting for this puzzle is that Fe(Te,Se) 
could be in a topological superconducting (SC) state. If it is 
the case, we can expect that nontrivial topology can integrate 
into the SC states in iron-based superconductors. 

Recently, some studies[7, 8 ] have revealed that the band 
structures can be tuned to have nontrivial topological proper¬ 
ties in monolayer Fe(Te,Se) and monolayer FeSe/STO. Fur¬ 
thermore, in electron-doped monolayer FeSe/STO, the exper¬ 
imental measurements have observed high temperature super¬ 
conductivity with Tc over 65K[9-16]. In analogy to the doped 
topological insulators, which are strongly believed to be topo¬ 
logical superconductors [4, 17-19], a natural question arises, 
can the electron-doped monolayer FeSe/STO be topological 
superc onductors ? 

In this paper, we propose that the electron-doped monolayer 
FeSe/STO could be an odd-parity topological superconductor 
in the spin-triplet orbital-singlet s-wave pairing channel. To 
show this exotic state, we first analyse the distinctive elec¬ 
tronic structure of monolayer FeSe/STO, and present a sys¬ 
tematic classification of the pairing symmetry in monolayer 
FeSe/STO from the lattice symmetric group. Second, we dis¬ 
cuss the topological properties of such odd-parity SC states, 
and extract the minimum effective models to capture the es¬ 
sential physics. Third, we calculate the phase diagram of SC 
states according to different scenarios of effective pairing in¬ 
teraction. Finally, we discuss the experimental signatures of 
the topological SC states. 



FIG. I: (a) The Se-Fe-Se trilayer structure. The black/green balls 
with deep and light filling label Fe/Se atoms. Here, the deep/light 
filling Se atoms are above/below the Fe plane. The red/black dashed 
squares label the one-Fe/two-Fe unit cells, (b) The Fermi surface 
of monolayer FeSe/STO is schematically illustrated. The red/blue 
electron pockets have odd/even orbital parity. The red/black dashed 
squares label the one-Fe/two-Fe Brillouin zone.The evolution of the 
band structure from (c) the free-standing monolayer FeSe to (e) 
monolayer FeSe/STO with tensile strain. The red/blue color labels 
the spectrum with odd/even orbital parity. 


II. PAIRING SYMMETRY CLASSIFICATIONS 


The lattice structure of monolayer FeSe is shown in Fig. 
1(a). The two-Fe unit cell includes two Se and two Fe la¬ 
beled by A and B. The space group PA/nmra governs the 
Se-Fe-Se trilayer structure, and belongs to a non-symmorphic 
group[20-23]. Indeed, there exists a n-glide plane described 
by the operator which involves a fractional trans¬ 

lation ( 5 ^) combining with the afo-plane mirror. Centered on 
an Fe atom (see Fig. 1(a)), eight point group operations, E, 
2 S' 4 , C 2 {z), C 2 {x), C 2 {y) and 2ad, form a D 2 d point group. 
Together with an inversion followed by fractional transla- 
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tions ( 55 ), i.e., they generate all the elements of 

P4:/nmm. The 16 operations do not form a point group. 
However, if the fractional translation (^ 5 ) is stripped off, the 
16 operations form a point group, which indeed is Dih- It 
is convenient to classify the pairing symmetry with the ir¬ 
reducible representation (IR) of £> 4 / 1 . For this purpose, one 
simple way is to recompose the Bloch wave functions in the 
one-Fe Brillouin zone (BZ). 

The glide plane symmetry divides the five d or¬ 

bitals into two groups {dxz,dyz) and {dxy,dx 2 -y 2 ,dz 2 ), and 
each group is recomposed to be the eigen-states of the glide 
plane operation with the definite orbital parities. The tight- 
binding Hamiltonian can also be decomposed into two parts 
with inverse orbital parities, which allow us to transfer the 
two-Fe unit cell picture into one-Fe unit cell picture[20-22]. 
In momentum space, the tight-binding Hamiltonian in one-Fe 
unit cell picture can be written as 

ijo = E +Y, rJimeikWaO^). 

k.cr k,cr 

( 1 ) 

Here, the first/second term has odd/even orbital parity un¬ 
der the glide plane operation. 1 p^{k) = [dxz,aO^), dyz,<jO^), 
dx2-y2^cr0<-), dxy,aOi-), 42 _^(k)]'^ withd™,cr(k) denoting the 
electron annihilation operator at the mth orbital with momen¬ 
tum k and spin a. = '^^(k + Q) and Ae(k) = 

^o(k -f Q) with Q = (tTjTt) (see Appendix A for details). 
The energy spectra from Eq.(l) are shown in Fig. 1, in which 
Fig. 1(e) is consistent with observations of the angle-resolved 
photoemission spectroscopy (ARPES)[10, 11], and the chem¬ 
ical potential is set to satisfy that 10 % electrons is doped per 
Ee clarified by experiments[10-12]. The fundamental dif¬ 
ference between Eig. 1(c) and (f) is referred to the band- 
renormalization effect induced by the strain from the STO 
substrate, which strongly modulates the hopping parameters 
between the {dxz, dyz,dxy) orbitals and switches the positions 
of two doubly degenerate points Mi and M 3 at Mx high sym¬ 
metric point, where the Mi point mainly has {dxz,dyz) or¬ 
bital weight and the M 3 point mainly has dxy orbital weight. 
This picture is the most natural and simplest to account for the 
distinctive electronic structure of monolayer EeSe/STO com¬ 
pared to other scenarios [24—26]. 

The SC order parameters should follow the IRs of the sym¬ 
metry group of the system. It is safe to use Dih to do so in the 
picture of one-Ee unit cell according to our aforementioned ar¬ 
guments. There exist two kinds of symmetry-allowed Cooper 
pairs, i.e., (k, —k) and (k, —k + Q) pairing channels. Previ¬ 
ously, the (k, —k-f Q) pairing channels are proposed to coex¬ 
ist with (k, —k) pairing channels to explain the nodeless and 
sign-change gap structures in iron-based superconductors [ 20 , 
21]. The price for coexistence of both kinds of pairings is 
that the orbital parities are mixed and the spatial inversion 
symmetry is broken. Here we focus on an SC state with 
only one IRs in (k, — k) pairing channel and leave to discuss 
the irrelevant (k, —k -|- Q) pairing channel in Appendix B. 
Moreover, we only need to consider the pairings between the 
three t 2 g orbitals as the orbital weight for Eg orbitals are ne- 
glectable on the Eermi surfaces[27]. Define the Nambu basis. 


TABLE I: The IRs of all the possible onsite superconducting pairing 
in (k, —k) channels. Here, 771/4 = =p|(Ao + 2'\/3A8) and 772/3 = 
|(=FAo ± %/ 3 A 8 =F 3A3/1). 
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T'(k) = [{dt(k)},{d;(k)},{d|(-k)},{- 4 (-k)}]'^ with 
{d^(k)} = {dxz,a{^),dyz,a{^),dxy,c,{\^)}. The pairing term 
in the Bogoliubov-de Gennes (BdG) Hamiltonian can be ex¬ 
pressed as 

= (^2) 
k 

Here, Tx is one Pauli matrix in Nambu space, and A(k) is 
a 6 X 6 matrix. Our purpose is to identify the exact form 
of A(k). Eor convenience, we utilize four Pauli matrices 
(sq, Sx,Sy, Sz) to span spin space and nine Gell-Mann matri¬ 
ces (Aq, ..., Ag) (see Appendix B for definitions of Gell-Mann 
matrices) to span orbital space. In such a way, A(k) can be 
decomposed into the product of the Pauli matrices and Gell- 
Mann matrices, i.e. A(k) = /(k)smA„ in which /(k) is the 
pairing form factor. We summarize all the possibilities of the 
(k, —k) on-site pairing channels according to the IRs of 
in Table I and non-on-site pairing channels up to the next- 
nearest neighbor in Table II. In both Tables I and II, the spin- 
singletZ-triplet pairing channels are listed in the first/second 
parts. 


III. TOPOLOGICAL SUPERCONDUCTING STATES 

To evaluate the pairing channels that could support the 
topological SC states, we first impose the nodeless gap struc¬ 
tures restrictions to the pairing channels in Tables I and II 
according to ARPES and STM/S experimental results[9-l 1], 
i.e., Eu \ Eu^ and E^ ^ in Table I and with 

/ 4 (fc)soAo, , a\2 and E^'^ in Table II. Second, we focus 
on the odd-parity pairing channels based on the proposals that 
odd parity pairings usually support the topological SC states in 
doped topological insulators [4]. Einally, we consider the SC 
states with the C 4 rotation symmetry verified by both exper¬ 
imental observations[10-13] and our calculations in Section 
IV. This constraint forces the time-reversal (TR) symmetry to 
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TABLE II: The IRs of all the possible nearest and next nearest neigh¬ 
bor superconducting pairing in (k, —k) channels. Here, fi/ 2 {k) = 
coskj: ± cosfcj,; f4{k) = cos k^o cos ky-, [fsik^), fs^ky)] = 

[sin fca;, sin fc^j; / 5 (fc) = sin sin ky. 

(k,-k) : A(k) iR“ 

fl/4{k)soXo/s, f5{k)soXl, f3{k:c)soX5 -I- f3{ky)soX7 A\J 

/2(A:)soAo/8,/3(A:a:)soA5 - f3iky)soX7 

f2{k)soXl, f3(ky)soX5 — f3{kx)soX7 A2g 

/5(fe)soAo/8, /l/4(fe)soAl, /3(fcy)soA5 + f3{kx)soX7 B 2 g 

^/i/4(^)^zA2, ^ ^ ^ \f3{kx')Sz/x/yXi 3“ fsiky^^z/y/xX&\ A\g 
if2{k^SzX2^i ^ ^ [f3{kx')Sz/x/yX4 ’^fsi^ky'jSz/y/xXo] ^Ig 

^ ^ ^ [/a )5z/x/y A4 ^3 (^3: )^2/y/x As] ^2g 

^/s (^)'5z A2 , I ^ ^ 1/3 (^i/)'53/a;/y A4 -|- /3(^x)Sz;/y/3;A6] .^2g 

‘ifl/2/4/5{k){Sx, Sy)X2 Eg 

f3(^kx')Sx/yXo i f3{ky^Sy/xXQ 

[/3(fcx),/3fe)]gzAo_^ 



be broken spontaneously for some states. With all above 
constraints and turn to the monolayer FeSe/STO, four possible 
odd-parity pairing states survive: (1) E^\ a doubly degener¬ 
ate TR breaking state with Ai(k) = AoS 2 (A 5 ± 1 X 7 ), (2) 

/n't 

Eu , a TR invariant state with A 2 (k) = Aoz(sa;A 5 + SyXr), 
(Note that E^ ^ is equivalent to E^^), (3) E^\ a doubly 
degenerate TR breaking state with A 3 (k) = Ao[/ 3 (fca;) ± 
if 3 {ky)]szXo, (4) a TR invariant state with A 4 (k) = 
Aq [/s (fca;) Sa; Aq +/3 (fcy) Sy Ao], (Note that all four components 

in h{kx)sx/yXo ± h{ky)Sy/xXQ} are equivalent.). 

Through the bulk-boundary correspondence, we demonstrate 
that all these four kinds of odd-parity pairing channels sup¬ 
port topological SC states in monolayer FeSe/STO. The BdG 
Hamiltonian describing the SC states can be obtained by com¬ 
bining the tight-binding Hamiltonian Hq in Eq. (1) and pair¬ 
ing term Hp in Eq. (2), i.e., 

HBdG = Ho + Hp. (3) 

Note that HBdG in Eq. (3) includes both odd-orbital-parity 
and even-orbital-parity parts. The edge spectra from the odd- 
orbital-parity parts of ElBdG with Ai(k)...A 4 (k) are pre¬ 
sented in Eig. 2. The even-orbital-parity parts of ElBdG give 
the same spectra if ky is translated to ky + tt (see Eig. 1(b) 
for comparison). The edge spectra in Eig. 2 explicitly sup¬ 
port the Andreev bound states which are the identifications 
of topological superconductors. Besides, the bulk proper¬ 
ties of topological superconductors are usually characterized 
by some topological numbers. Here, the pairing channels 
with Ai(k) and A 3 (k) break the TR symmetry, and the 
Chern number[28] can be introduced to characterize such two 

states, i.e., C = / dk{VkUn(k)\ x |VA:Wn(k)). 

JBZ 

The calculations show that both odd-orbital-parity and even- 
orbital-parity parts give the Chern numbers C° = C® = 4 in 
the one-Ee BZ for Ai(k) and A 3 (k) pairing channels. Thus, 
such two pairing channels are characterized by the total Chern 
number C =^(C°+ C®) = 4 in the two-Ee BZ. The Chern 
number C = 4 is equal to the number of edge Andreev bound 


FIG. 2: The edge spectra of odd-orbital-parity BdG Hamiltonian with 
Ai(A:), A. 2 {k), A 3 {k) and A 4 (fc) in (a), (b), (d), and (e). In the 
presence of the orbital-parity-broken perturbation, i.e., the staggered 
potential of Fe sublattices, the edge spectra of BdG Hamiltonian with 
A 2 {k) and A 4 {k) are shown in (c) and (f). Here, the system has pe¬ 
riodic boundary condition along the y direction and open boundary 
condition along the x direction with 51 one-Fe unit cell lengths. The 
red/blue colors label the edge states localizing at the opposite bound¬ 
aries, and the dashed/solid lines label the edge states with up/down 
spin directions. Note that the degenerate edge states on the same 
edge are artificially split as a guide for the eye. 


states shown in Fig. 2 (a) and (d). For the TR invariant A 2 (k) 
and A 4 (k) pairing channels, the total Chern numbers are zero. 
However, the spin Chern numbers[29, 30] can be introduced 
to characterize the bulk topological properties of SC states 
in A 2 (k) or A 4 (k) pairing channels. Namely, = 1, 

= — 1 in the two-Fe BZ. Correspondingly, two Z 2 topo¬ 
logical numbers[31] with opposite orbital parities defined by 
^o/e _ = 1 characterize the bulk topological 

properties for SC states in A 2 (k) or A 4 (k) pairing channels. 

Having confirmed that the topological SC states emerge in 
the nodeless odd-parity pairing channels, we notice that the 
edge spectra shown in Fig. 2 (a) and (b) and the edge spectra 
shown in Fig. 2 (d) and (e) are very different. Therefore, it is 
necessary to extract the minimum effective models to clarify 
the essential physics hidden behind. First, we are aware of 
the A 3 / 4 (k) pairing channels being in the intra-orbital spin- 
triplet p-wave pairing channels. Thus, the orbital degree of 
freedom is inessential, and the minimum effective Hamilto¬ 
nian can be reduced into the single band space, which is the 
same Hamiltonian to describe the well-known p±ip topolog¬ 
ical superconductors/superfluids [1, 32, 33], and the nontrivial 
topology is referred to the p±ip pairing terms. Therefore, we 
omit our discussions for these “trivial” topological SC states. 

For Ai(k) and A 2 (k), which are the inter-orbital spin- 
triplet s-wave pairing channels, the three t 2 g orbitals are in¬ 
volved and entangled with each other not only in the bands 
around the Fermi surface shown in Fig. 3 (a), but in the pair- 
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ing terms shown in Fig. 3(d). Note that we should have three 
bands when we consider three t2g orbitals. It indicates that the 
third band mainly with the dxz and dy^ weight has to strongly 
couple with two Cg orbitals and be gaped and pushed away 
from the Fermi level. In order to describe the two bands in ex¬ 
act three orbital basis, we adopt the angular momentum repre¬ 
sentation characterized by the azimuthal and magnetic quan¬ 
tum numbers I and m. The new electron creation operators are 
4 ™= 2 .±i).<t( 1 ^) = T-^[dlz,A'^)±idl,^^(k)], then we have 
A\{k) - 

Aj(k) ~ .^(k)4y,t(-k) -f 

Now, we can only exploit the operators involving in Ai/ 2 (k) 
to construct the basis to write the minimum effective Hamil¬ 
tonian, and this approximation is equivalent to treating with 
dxz and dyz orbitals with equal weights. In the effec¬ 
tive basis, T'i/ 2 (k) = [{■!/'i/ 2 t(k)}, {■i/'i/ 24 ,(k)}]'^ with 

{V'i/2.a(k)} = {d(2,i/-(-i)-),<T(k), dxy,a(k), 4y.ff/a(-k), 

Here, kis measured from the M point, a = —a, and (—1)'^ = 
1/—1 for spin the orbital parity index is omitted for sim¬ 
plicity. =r,[4^/^)(k)+n=.d!'/^^(kK] + 

TxAo, 'H^^^(k) = n[^\k) und'U'iHk) = ■Hf^*(-k). The 
three Pauli matrices (Ti/ 2/3 are introduced to span the effective 

two-band space, = ^i(k)+^ 2 (k) _ di^^^\k) = 

TAky, = -Akx and 

HA)(k) breaks TR symmetry, because only m = 1 is in¬ 
volved. (k) is TR invariant, and characterized by the 
T~^H^‘A{k)T = H'A) (—k), where the TR symmetry oper¬ 
ator is T = iSyToaolC with /C the complex conjugated op¬ 
erator. The dispersions £ 1/2 (k) with definite orbital parity 
can be read out from Fig. 1 (e) and Fig. 3 (b). Around My 
point, we have e\i^{k) = exji ~ /^ + + /5i/24 

£°y 2 (k) = ei /2 - -f / 3 i/ 2 fcs + ai/24' 'The signs of a//3 
are crucial to determine the properties of the topological SC 
states. In Fig. 3 (b) and (c), we schematically illustrate the 
evolution of the £°y 2 (k) under the couplings induced by glide 
plane around My point, and we can find ei < 62 , ai < 0, 
Pi > 0, Qf 2 > 0, P 2 < 0. The effective mass measuring 
the energy gap Em^, — Emi shown in Fig. 1 (e) or Fig. 3 
(c) is m = '^^ 2 '^^ > 0. The finite electron-doped condition 
/i^ -f > to^[ 34] always supports topological SC states 
for where the chemical potential /i is measured 

from the middle of gap. The remarkable feature of the edge 
spectra in Fig. 2 (a) and (b) is that the edge Andreev bound 
states have a twist (three times of crossings) around ky = tt 
and only one crossing around ky = 0. This difference can 
be understood with the “orbital mirror helicity” from the mir¬ 
ror operator in C2{x/y) acting on three t2g orbitals in analogy 
to the “spin mirror helicity” proposed in Ref [34]. The con¬ 
servation of mirror helicity force the non-/twisted feature of 
the edge Andreev edge states under the non-/band-inversion 
conditions between £^'^°(k) and £ 2 ^°(k) along x direction. 



FIG. 3: (a) The weight of three t 2 g along the Fermi surface around 
My with odd-orhital parity, (h) and (c) The effective band disper¬ 
sions without/with inter-orbital coupling from glide plane, (d) Three 
competitive pairing channels with (/? = ^in weak-coupling limit. 


sgn[(e 2 -ei)(a 2 -ai)] > 0 /sgn[(e 2 -ei)(/ 32 -/ 3 i)] < 0(Note 
that e ®^2 (My -F k) = £“^2 (Ma; -F k)). We are aware of the im¬ 
portance of the non-symmorphic lattice symmetry which not 
only induces the orbital-momentum locking k x cr • z through 
glide plane, but protects the exotic behaviors of the edge An¬ 
dreev bound states. We can verify this point through intro¬ 
ducing the staggered onsite potential, which mixes the orbital 
parities, breaks the non-symmorphic lattice symmetry and de¬ 
stroys the twist feature of edge spectra. The results are shown 
in Fig. 2 (c) and (f). However, the bulk topological properties 
are robust against such perturbations. 


IV. THE EFFECTIVE PAIRING INTERACTIONS 

Although the high temperature interfacial superconductiv¬ 
ity in monolayer FeSe/STO seems to have been established 
beyond doubt, the mechanism for superconductivity is still 
an open question[35], and the unique features of monolayer 
FeSe/STO further pose a higher barrier to block our under¬ 
standing of the superconductivity from some standard theo¬ 
ries. For example, the monolayer FeSe/STO is strictly two- 
dimensional and has no hole pockets at the BZ center while 
its three-dimensional counterpart bulk FeSe resembles iron- 
pnictide with hole pockets. The Fermi surface of mono- 
layer FeSe/STO is similar to that of A 3 ;Fe 2 _ySe 2 (A=K, Cs, 
Rb), except that the small electron pocket around (0,0, tt) 
in Aa;Fe 2 -j,Se 2 is absent here. In weak coupling limit, 
the spin-fluctuation-exchange theory predicts that the {Big; 
f 2 {k)stjXo} pairing channel is dominant in A 3 ;Fe 2 -gSe 2 and 
the gap structure has nodes along the kz direction[36, 37]. 
However, the ARPES measurements reported isotropic full 
gaps without nodes on all pockets in A3;Fe2_ySe2[38, 39]. In 
the strong coupling limit, the phenomenological t-J model 
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predicts that the {Aig\ / 4 (fc)soAo} pairing channel is domi¬ 
nant in A 3 ;Fe 2 _ySe 2 and the gaps have same sign for all the 
pockets [40]. However, the inelastic neutron scattering mea¬ 
surements on A^Fe 2 _ySe 2 reported a resonance with wave 
vector Qc = (7r,7r/2) in the superconducting state[41], 
which indicated that there existed a sign-change between 
the fermi surfaces connected by Qc. These contradictions 
strongly question the standard theories. On the other hand, 
the studies of some confirmed systems with interfacial su¬ 
perconductivity including bilayer lanthanum cuprate[42] and 
LaAlOa/SrTiOa heterostructure [43] could provide us some 
useful insights to understand the superconductivity in mono- 
layer FeSe/STO. The studies of the aforementioned systems 
indicate that surface phonon plays a key role to drive the 
superconductivity [44]. Recent ARPES experiment observed 
the band replication, which was attributed to strong coupling 
between the cross phonon and electrons [15], and the cooper¬ 
ation between the cross phonon mode and spin fluctuation is 
argued to be the origin to enhance Tc in monolayer FeSe/STO. 
Therefore, it is still possible that the superconductivity in 
monolayer FeSe/STO is driven by the electron-phonon cou¬ 
pling, and the surface phonon-mediated SC mechanism in 
monolayer FeSe/STO has been proposed in Ref. [45]. Here, 
without loss of generality, we consider several possibilities of 
the effective interactions that can drive superconductivity in 
different pairing channels and focus on the parameter regime 
missed previously. 

We first assume the multi-orbital Hubbard interactions as 
pairing driver. 


i,l 

Jh ^ ^ (2S2/ ■ “h 

z,Z>Z' 

4- J ^ (5) 


Here, U, V, Jh, J' are the intra-orbital, inter-orbital, Hund’s 
coupling and pairing hopping term. I, I' G {xz,yz,xy), 
and Sii = ^dh^Saa'diia- The spin rotation symmetry re¬ 
quires U = V + 2Jh, and Jh = J' at the atomic level. 
Since the predictions from the weak-coupling theory [36, 37] 
about iTo+ Hint were not consistent with the experimen¬ 
tal reports[38, 39], the strongly correlative picture with quite 
large Jh is possible and the strongly correlative effects in 
iron chalcogenides have been reported by recent ARPES 
experiments [46]. Deflne the pairing operators. 


k k 


^ -f (-k)], 

crcr' 

Kwik) = - dl,n{k)dtn'i-k)]. 

gg' 

( 6 ) 

The interaction Hamiltonian has the form, 

= [/ ^ kluKu + JhJ2 

I liii' 

+ {V-Jh)J 2 + iv + Jh) 

ll'a 1=^1' 

(7) 

When the Hund’s coupling is strong enough, i.e., Jh > U /3, 
the third term of Eq.(7) can give rise to the instability in a spin- 
triplet channel[47, 48], which involving the {Aig: 182X2}, 
and IRs in Table I. The detailed discussions about 
these pairing channels are merged into the third kind of effec¬ 
tive interaction in the following. 

Another standard theory for the superconductivity is the 
phenomenological Heisenberg model in strong coupling 
limit, we consider the effectively frustrated Heisenberg 
interaction [49] as the pairing force, 

= E S*z-S,, + J2 E Szi-S,,. (8) 

Here, J1/2 are the nearest and next nearest neighbor mag- 

( 2 ) 

netic exchange couplings. A well-know result of H^J^ is that 
the magnetic ground state is checkerboard-antiferromagnetic 
when 2 J 2 < |Ji|, and collinear-antiferromagnetic when 

2 J 2 > |Ji|. However, no Eermi surface reconstruc¬ 

tion induced by spin density wave was observed in mono- 
layer EeSe/STO but in mutli-layer EeSe/STO in ARPES 
experiments [12]. The recent first-principles calculations pro¬ 
posed that the magnetic order was strongly frustrated in mono- 
layer EeSe/STO with 2 J 2 Ri |Ji|[50]. Another issue is 
the sign of Ji. If both Ji and J 2 are antiferromagnetic, the 
A 3 / 4 (k) pairing channels are ruled out, and the SC states fall 
into {Aig-. /4(fc)sQAo} induced by J 2 or {Big-. / 2 (A:)soAo} 
induced by Ji. If Ji is ferromagnetic and J 2 are antiferromag¬ 
netic, the A 3 / 4 (k) pairing channels are possible from sym¬ 
metry point, but these two odd-parity pairing channels have to 
compete with the {Aig-. f4{k)soXo} induced by J2. The win¬ 
ner is determined by the topology of the Eermi surface[51]. 

Eor the low electron-doped at O.le/Ee, the Eermi pockets lo¬ 
cating at M points are quite small. Therefore, the form fac¬ 
tor fiik) has large magnitude, and the SC states favor the 
{Aig-. / 4 (fc)soAo}. Ifthe electron-doped level can be tuned in 
monolayer EeSe/STO without suppressing the superconduc¬ 
tivity. We can expect that the SC states in over electron-doped 
samples would favor A 3 / 4 (k) pairing channels for ferromag¬ 
netic Ji, because the Eermi surface locates at the X points, 
where the form factors fsik^/y) have large magnitudes. We 
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note that such kind of pairing was discussed in underdoped 
cuprates [52]. 

From the aforementioned arguments about the possibly 
significant role of surface phonon, we consider the third 
kind of phenomenological interaction from phonon-mediated 
mechanism[45] to induce the interfacial SC instability in 
monolayer FeSe/STO, 


rr(3) _ 
^int 


^ 2 

1,1' ,cr,cr' k.k' 




k,l(T 


fd—k'Va'^k'^la- 


(9) 

Here, we assume (k, k') = —Vq for I = V, a' = a and 
(k, k') = —Vi for I > V. Note that the third term in 
Eq.(7) with Jh > C//3 can also be described by With 
the pairing operators shown in Eq. (6), takes the form. 


= -V, E KiAs,ii - E 


i>i' 




( 10 ) 


l>l'a 


Under the mean-field approximation, Ag/;/ = (A^ ;j,), 

= (A“;j,), the 1*® decoupled as follows, 

= -Vo E E 

I l>l' 

- Ul E + H.C. + heon. (11) 

Here hcon = J2iVo\^s,ii\'^ + P + 

Vi a I Ag i/' P- Now, we consider the odd-orbital-parity 
parts of the normal-state Hamiltonian. The mean-field Hamil¬ 
tonian takes the following form. 


Hmf = E + H,on, (12) 

k 

where T'(fc) has the same form 
shown in Eq.(2) except {(icr(k)} = 

{_dxz,a{}V) ■) dyz^a (k), dxy^cT (k), dx^ —y2 fj (k), dz2 (k) } 

now. Then HMF{k) = Ho{k)Tz + A{k)Tx, 

5 

= Ao{k) 0 Ao{k), and Hcon = J2 ^ o,mm (.k) + hcon- 

k,m—l 

Assume the HMrik) can be diagonalized with matrix 
Uk, i-e., Ul.HMF{k)Uk = Ek,i ® -E'fc. 2 ----E'fc. 20 - Then the 
mean-field self-consistent equations take the forms. 


Fk,n,fik,n,l' + 10 + U^^^l_^_^Uk,n,l' + 15\f{Ek,n) 

2 ’ 

— {^k,n,fik,n,l' + 15 + t^fc,r!,,i+5^fc."7' + lo]/(^fc,ra) 

2 ’ 

-'^{Uk,n,Fk,n,l' + 15 “ Uk,n,l+5^k,n,l' + w]f{Ek,n) 
2 

— [^k,n,filc,n,l' + 10 ~ U^^^ i^^Uk,n,l' + lo]fiEk,n) 

2 ’ 

10 

Ae = E E \U*k,n.J^fiEk,n). (13) 

k,n—l m—1 

Here, /(x) = — is the Eermi distribution function 

e'=B'-'' +1 

and Nf, is the electron number. In comparison with Ta¬ 
ble I and Eq.(ll), the relevant IR channels in Table I can 
be represented with (13). Eor example, soAo}= 

^ 2 ) 

■^0 (^ 5 , 3 : 2 :, 3 : 2 : 0 Ag^yz^yz 0 Ag^xy,xy ), {Eu : i{sxX5,SyX7)}= 

^i^la:z,xvSxX5, Al^^^^^SyXr). Likewise, other IR channels 
can be read out following same way. 

It is possible for A(/c) to take the form of linear combi¬ 
nations of several different IR channels, but some symme¬ 
tries have to be broken to pay the price for such coexis¬ 
tence. Eor example the inverse symmetry is broken for the 
SC states proposed in Refs. [20, 21]. Likewise, the TR sym¬ 
metry or lattice symmetry could also be broken when two 
different one-dimensional IRs or two components in a two- 
dimensional IR coexist. In order to gain some insight be¬ 
fore we perform the numerical calculations, we note that all 
the experiments reported the isotropic Eermi surface and gap 
structures without any resolvable distortions, and the mono- 
layer EeSe/STO was conformed to be the cleanest composi¬ 
tion with the simplest structure[10-12]. These features rule 
out the possibilities of some complex orders, such as nematic 
order found in bulk EeSe. Erom Table I, we can first elim¬ 
inate the possibilities of the {B 2 g'. sqAi}, {Aig- iSzX 2 }, 
{Eg-. i{sx,Sy)X 2 } and {Aig\ so^g} pairing channels, be¬ 
cause the leading vi\X.&!:-dxz-dyz hopping term is proportional 
to sin kx sin ky, which is nearly zero around the Eermi surface, 
and the {Aig: spAg} channel has nodes. Second, it is straight¬ 
forward to check that two components in {Eu- so(A 4 , Ag)} or 
{E^k_ S 2 ;(A 5 , At)} give two degenerate strip SC states with 
nodes. Thus, the TR-broken linear combination of two com¬ 
ponents is optimal to achieve the isotropic nodeless gap struc¬ 
ture and lower the energy. Note that the coexistence of these 
two two-dimensional IRs could rise the energy, because they 
follow different transformations under the lattice symmetric 
operations and suppress the gap amplitude. Einally, no ad¬ 
ditionally global symmetries can be broken for sqAo} 

(‘2'\ 

and {£’i . z(sa:A 5 , SyAr)} to coexist with each other and with 

{Eu- so(A 4 , Ae)} or {E^^: S 2 (A 5 , A 7 )} to avoid breaking the 
isotropic SC gap structure and achieving lower energy. There- 


20 


A..zi' = E 

k,n— 

20 

Amp = E 

k,n— 

20 

Kw = E 

k,n— 

20 

Amp “ E 


k,n—l 

20 
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FIG. 4: (a) The evolution of three components of SC order parame¬ 
ters in channel about Vq. (b) The evolution of components of 

SC order parameters in Eu, and E^'’ channels about Vi. (c) 
The evolution of the condensed energy in different SC states with 
relevant IRs about Vo and Vi. (d) The phase diagram is plotted in 
(Vo, Vi) plane with respect to the lowest energy. We set 51 xSlmesh 
of k, and the electron number to satisfy electron-doped O.le/Fe. The 
energy scale is measured with eV. 

fore, we find that these four IRs, i.e., {E^. So(A 4 , Ag)}, 

^^(As, At)}, soAqI and {E^^\ i[sx\b, Sj/At)} are in¬ 

dependent, and TR symmetry should be spontaneously broken 
in the hrst two IRs. It is straightforward to verify these argu¬ 
ments through the following numerical calculations. 

Now, we perform the numerical calculations to evaluate 
which pairing channel governs the ground state of the sys¬ 
tem for different Vq and Vi. The ground state energy of 
Eq.(12) is Gs{T) = -fesTlnTre-^^"^, and Gs{T ~ 0) = 
10 

{Hcon ~ ^ S \^k,n\) zcro temperature. For simplic- 

k,n—l 

ity, we can evaluate the ground state through the minimum 
of the condensed energy density dehned as fg = hcon — 

10 5 

f - 4^7 E / d^kjE^^J for given electron 

n—1 n—1 

number, where E^ ^ are the energy spectra of normal state. 
Solve the self-consistent equations (12) and (13) for param¬ 
eters (Vo,Vi) with respect to the minimum of fg, we show 
the evolution of SC order parameters and condensed energy 
about (VqjVi) in Fig. 4, and we hnd topologically triv¬ 
ial soAqI channel and topologically nontrivial {Ed- 

i{sxX 5 ,SyXr)} are dominant in relevant regime of (Vq, Vi) 
parameter plane. 


V. DISCUSSION AND SUMMARY 

If the superconductivity in monolayer FeSe/STO is driven 

(3) 

by the effective interaction in Fq. (10), the observed 


isotropic and nodeless s-wave gap structures select both topo¬ 
logically trivial (sqAo) and nontrivial i7}^\A2(k)) as 
possible candidates. The essential difference lies in that the 
former one has even-parity and spin-singlet pairing while the 
latter one has odd-parity and spin-triplet pairing. Therefore, it 
is unambiguous to adopt the experiments which can directly 
distinguish the spin states and parities to pin down the possi¬ 
ble candidate. Particularly, temperature dependence of the nu¬ 
clear magnetic relaxation(NMR) rate can be utilized to distin¬ 
guish the two different pairings. The well-known result is that 
the NMR rate has a Hebel-Slichter peak at the SC transition 
temperature for the even-parity and spin-singlet s-wave SC 
state[53]. However, the Hebel-Slichter peak could disappear 
with the anti-peak behavior due to the unique spin, orbital, 
and momentum locking effect in topological SC states with 
odd-parity as shown in Ref. [54]. The parity of the Cooper 
pair is characterized by the inverse operator {*|^^}. It indi¬ 
cates the odd-parity pairing has a sign change or phase shift 
of TT between the top Se and and bottom Se layers along the 
c-axis compared with the even-parity pairing. Thus, the stan¬ 
dard magnetic-flux modulation of dc SC quantum interference 
devices (SQUIDS) measurements [4, 55, 56] provide another 
scheme to distinguish the odd- and even-parity pairings. On 
the other hand, some transport measurements can also be ap¬ 
plied to detect the topological superconductors, such as the 
thermal Hall conductivity [57, 58]. The challenge for such 
measurements is that the FeSe is very air sensitive, and the 
experimental measurements should be performed under the 
ultra-high vacuum condition. 

In the aforementioned discussions about the SC pairings, 
we assume that the glide-plane symmetry is not broken. Ac¬ 
tually, there exist some possible effects to break the glide- 
plane symmetry. For example, the atomic spin-orbital cou¬ 
pling could be have non-neglectable effect in iron chalco- 
genides. It is explicit that the inter-orbital spin-orbital cou¬ 
pling can mix the bands with inverse orbital parities, and in¬ 
duce the inter-orbital SC pairing in (k, —k + Q) channels. 
However, the weight of intei-dxz-dxy spin-orbital coupling is 
proportional to Ago ~ 0.05eU[37], while the intei-dxz-dxy or¬ 
bital hopping term with definite orbital parity is proportional 
to \2itlfsmkp\ ~ 0.3eU at the Fermi surface. We can es¬ 
timate that the ratio between the SC pairing order parameter 
in (k, —k + Q) channel and that in (k, —k) channel should 
be ^ 0.025. Thus, the atomic spin-orbital coupling plays a 
neglectable role in SC states, and the Ref. [37] verified this 
point in Aa;Fe 2 _ySe 2 . Other issues, such as the coupling 
between the monolayer FeSe and substrate STO, could also 
break the glide-plane symmetry. Such couplings are tunable 
and strongly affected by the fabrication process and the sub¬ 
strate materials[13, 59]. Here, we consider the case that the 
strength of coupling between the monolayer FeSe and sub¬ 
strate is weak in comparison with the relevant hopping ampli¬ 
tude. 

Compared with the general topological materials, in which 
the extended s and p orbitals are the bricks to build low-energy 
electronic structures, and the spin-orbital coupling plays an 
essential role in inducing the strong linear couplings, the lin¬ 
ear couplings in monolayer FeSe/STO is attributed to effective 
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couplings between 3d orbitals induced by d-p hybridizations 
from the unique non-symmorphic lattice structures. Such fea¬ 
tures provide us an alternative route to search for the new 
topological materials in strongly-correlative electronic sys¬ 
tems. 

In conclusion, we propose that the monolayer FeSe/STO 
could support the odd-parity topological SC states with the 
nodeless s-wave gap structures. In contrast with other topo¬ 
logical superconductors [2, 4]in which the spin-orbital cou¬ 
pling plays a key role, such topological SC states have strong 
relations with the unique non-symmorphic lattice symmetry 
which induces the orbital-momentum locking. Furthermore, 
we calculate phase diagram and suggest some experimental 
schemes to identify such uniquely nontrivial topological SC 
states. 
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Appendix A: The tight-binding Hamiltonian from symmetry 
analyses 

In this section, we discuss the properties of the tight- 
binding Hamiltonian from the symmetric point. The trilayer 
structure of the monolayer FeSe is shown in Fig. 1 (see main 
text). We focus on the three space group operations includ¬ 
ing glide plane symmetry operator, = {m^lro} with rg = 
( 5 ^) and two reflection symmetry operations, = {mx\ro} 
and gx' = {mj^/lOO}. Besides, the lattice has inverse symme¬ 
try denoted by the operator, gi = {i|ro}. According to the 
LDA calculation, we can only focus on Fe atoms, the Bloch 
wave functions for the 3d orbitals of Fe are defined as 

\o^V, k') = </'«(!■' - r^r,)- (Al) 

^ n 

Here, -f with lattice vector and the position 


with 


(A4) 


(k') = [4.xz(k'), d^y^ik'), 4,xy(k'), 

In the basis fl'(fc'), the corresponding transformation matrices 
for the three operations g^ have the following forms. 


Uigx') = 

rux' 0 

0 nix' 

Uigz) = 

0 

g-i(mjk')-ro 

Uigx) = 

0 

g-i(m^k')-ro 


g-i(m,k')ro 


niz 


0 

=k': 

0 


g-i(m,k')ro 


(A5) 


Where 


■ 0 

1 

0 

0 

O' 


- -1 

0 

0 

0 

O' 

1 

0 

0 

0 

0 


0 

-1 

0 

0 

0 

0 

0 

-1 

0 

0 

= 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

1 

0 

. 0 

0 

0 

0 

1. 


0 

0 

0 

0 

1. 


-100 0 0 

0 10 0 0 

mx= 0 0 10 0 

0 00-10 

0 0 0 0 1 


The symmetry of the Hamiltonian requires 

Hoik') = Uik')HoiUk')U\k') 

Define 


iJo(k') = 


HAik') HABik') 
HBAik') Heik') 


We can get 

HA/bO^x' ; ) — ‘^x' diA/B ( ^x' ; )'^x 

dd^AB i^x' 5 ^y' ) — dd ab ( d^x' i )’^x' 


(A6) 


(A7) 


(A 8 ) 


(A9) 


r'y of Fe atom r] = A, B, and (pa denotes the d orbital basis 
function {a = xz, yz, — y'^,xy, z^). The symmetry oper¬ 
ators acting on the basis function \ar], k') have the following 
properties, 

gx'\ari,k') = '^mx',a/3\l3r],mx’k') 

/3 

gz\ar],k') = ^ Wzk') 

/3 

gx\ar], k') = ^ '^'’"'mx^agl^rj, mxk'). (A2) 

/3 

The relevant tight-binding (TB) Hamiltonian can be expressed 


HAikx',ky') = m^HBikx',ky')m^ 


HaB 1 ky' ) 

= m^HBAikx',ky')m:, 

(AlO) 

^Asi^x'■> ky') = 

nixHBi-ky', -kx')mx 
nixHBAi-ky',-kx')mx. 

(All) 


Moreover, since \ar], k' -f G') = '’'i \ag, k') 

HA/Bi^' + G')=Ha/b{^') 

HABik' + G') = e'^'doHABik'). (A12) 

y'o = r'^ — = ( 5 , Considering the operator 5 ^, we can 

And in the entire BZ 


i?o = ^'f'^(k')i?(k')'I^(k'), 

k' 


(A3) 


0 


'HAik') HABik')]' 

niz 0 

5 

HBAik') HBik') \ 


(A13) 


as 
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We have 


0 rrix 

V = 

’ - 

dsxS 

0 

ruz 0 

- 

0 

hxb _ 

V-- 

1 

■ A 

A 


"71 

B 

-B 

1 


(A14) 


(A15) 


with A = l5x5,J3 = —rriz. It is straightforward to check that 
i/o(kO be block diagonalized, i.e.. 


The non-zero terms in A{k) are listed as follows, 

All/22 (*) = ei + cosfca: + COS fcy + cos fcx cos fcj, 

+ COS 2k^ + ‘^tfy/xx COS 2ky 

+ '^^lly/yyx COS 2fc^ COS ky + ^t]^yy/^^y COS kx COS 2ky 
+ '^t^Jxyy cos 2kx COS 2ky, 


V^Ho{k')V = iTii(k') 0 Tf22(k'), (A16) 

with iTii(k') = HAik') - iT^B(k')m 2 and i? 22 (k') = 
-ffyi(k') + HAB0^)'mz- From Eq.(A12), we can get 
^a/b{.^x* a ‘27TnxA ^x' 27r7Zy/) — HA/Bi.^x' A 

2 TTnx', kx' + 2 Trny') and HAB{kx' + ^nux', kx' + 2 'Kny') = 
^i{2-xn^,^+2Trnyi\)]^^B(kxi + 2'Knx',kx' + 2TTny'). When 
{nx',ny') = ( 0 , 1 ), iTii(k') = iTA(k') - HABik')m:, 
and H 22 {k') = HA{k' + Q') - Hab{^' + with 

Q' = (0, 27r). Furthermore, the momentum defined in the 
one-Fe BZ is kx = {kx' + ky')/2, ky = {—kx' + ky')j2 and 

Q = (7r,7r). 

Under the basis, T'l(k) = ['!/)l(k),'!/;l(k + Q)] 

with V'^(k) = [dl^{k),dl^{k),dl 2 _y 2 {k),dly{k),dl^{k)], 
di{k) = ^[dA,i(k') + (iB/(k')] and di{k + Q) = 

-^[(iA,i(k') - dB.i{k')] for I = xz,yz, di{k) = 
■^[dAji'ii-') - dB,i{^')] and d;(k + Q) = -^[dA,i{k') + 
dB.i{k')] for / = xy, x^ — y^, the TB Hamiltonian in the 
one-Fe BZ takes the following form 

iJo = ^4'^(k)iTo(k)4'(k). (A17) 

k 

Then, 

iTo(k) =iTo(k)0iJe(k). (A18) 

Here, He{k) = Ho{k + Q). 

The TB Hamiltonian in one-Fe BZ Eq. (A18) have block- 
diagonal forms, and each block has definitive orbital parity 
with respect to the glide plane symmetry. Besides, the in¬ 
version symmetry = {i|ro} indicates that the inversion 
center of monolayer EeSe is at the midpoint of Ee-Ee link. 
Thus, we can find that dxz/yz{'^)^dxy/x^-y^/z^{'^) are in¬ 
version even/odd, and dxz/yz{^ + Q,)/dxy/x^-y'^/z'^{^ + Q) 
are inversion odd/even. In other words, dxz/yz orbitals and 
dxy/x^-y^/z^ orbitals have opposite parities in the subspace 
with definitive orbital parity. The TB Hamiltonian in the one- 


Ee BZ is 

■ All A 12 

Ai3 

Ai4 

Ai5 


A 22 

A 23 

A 24 

A 25 

Ho{k) = 


A 33 

A 34 

A 35 




A 44 

A 45 

A 55 . 


A 33 (fc) = €3 0 2 f^^(cOS kx 0 cos ky) 0 dt^xy CO® kx COS ky, 


Ai4{k) =£4 0 2 f^"‘(cos kx 0 cos ky) 0 At'^y cos kx cos ky 
0 4 f^^j^(cos 2kx cos ky 0 cos kx cos 2ky) 

0 cos 2kx cos 2ky, 


J^55{k) = £ 5 , 


Ai 2 (fc) = — 4 f;).|sinfca; sinfcy, 


Ai 3/23 {k) = sin ky/^, ± 4itll sin ky/x cos , 


Ai4/24(fc) = -‘^itl'^sinkx/y + sin kx/y cos ky/x, 


Ai5/25(^) Anky/x sin ky/x cos kx/y. 


Ass{k) = 2t^{coskx — cosky). 


Ai5{k) = - sin kx sin ky. 

The onsite orbital energy is £i = £2 = 0.02. £3 = 
—0.539. £4 = 0.014 .£5 = —0.581, and the hopping pa¬ 
rameters for the free-standing monolayer EeSe are listed as 
follows[60], = -0.08/ - 0.311, = 0.232, = 

0.009/ - 0.045, = -0.016/0.019, t^lyy = 0.02, 

tf = 0.412, tfy = -0.066, = 0.063, ^4 = 0 .O 86 , 

4% = -0.017, ttiyy = -0.028, 41 = 0.099, = o.3, 

41 = -0.089, = 0.305, 41 = -0.056, 4^ = -0.18, 

41 = 0.146, ff = 0.338, tfy = -0.109. The renormal¬ 
ized parameters corresponding to Eig. 1 (d) in main text are 
tfy = 0.066, tf = 0.405, tf = -0.12. The renormal¬ 
ized parameters corresponding to Eig. 1 (e) in main text are 
44 ^ 0.076, tf = 0.183, tf = 0.405, tf = -0.311, 
41 = 0.19. 
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TABLE III: The IRs of all the possible onsite superconducting pair¬ 
ing in (k, — k + Q)channels 


(k,-k + Q) : 

A'(k) C 2 {z) 

C2ix) 

<Xd 

WU}' 

IR 

sqAo 

1 

1 

1 

-1 

Alu 

sqAs 

1 

1 

1 

-1 

Alu 

sqAi 

1 

-1 

1 

-1 

B2u 

so(A 4, Ae) 

(-1,-1) 

(1,-1) 

So(A6, A 4 ) (1,1) 

Ey 

iazX2 

1 

1 

1 

-1 

-Alu 

aziXs, X 7 ) 

(-1,-1) 

(-1,1) 

-az{X7,X3) (1,1) 

Ey 

i(^Sxj Sy)A2 

(-1,-1) 

(-1,1) 

i{ay,ax)X2 (-L-1) 

Eu 


(1,1) 

(1,1) 

-i{ayX7, ; 

SxAs) (1,1) 

Ey 

'^(SyAS; S3:;A7) 

(1,1) 

(-1,-1) 

-i{SxX7, , 

SyX5) (1,1) 

Ey 


TABLE IV: The IRs of all the possible non-onsite superconducting 


pairing in (k, —k + Q) channels 

(k,-k-fQ) : A'(k) IR^ 

/4,fcSoAo/8, /s.fcSoAi, /3,fc3,soA5 -I- /s.fcySoAr Aiu 

/2,feSoAo/85/3,fca,SoA5 — fi,kySo\7 -Blu 

/2,feSoAi,/s^fej^soAs —/s^fe^soAr A2u 

/5,fcSoAo/8, /l/4,fc'SoAl, f3,kySoX5 + f3,k^SoX7 _ B 2 u 

^fl/A,k^z^2^'i \f3,kx^z/x/y^4A fs^ky^z/y/x^o] Alu 

^y2,fe^zA2,i ^ ^ \f3,kx^z/x/y^A f3,ky ^ z / y / x^d] 

i^^°^°lf3,kySx/x/y^4 — f3,kxSz/y/xK] ^2u 

if5,kSzX2,i^^°^°[f3,kySx/x/y^4 + f3,kxSz/y/x^6] B2u 

ifl/2/4/3,k{ax, ay)\2 _ Eu 


The monolayer FeSe has inversion symmetry, thus every 
IR in Table I of main text should have a counterpart with 
a inverse parity. In other words, (k, —k -h Q) pairing 
channels should be possible from the symmetry point. For the 
(k, —k+Q) pairing. We define the the Nambu basis, T''(k) = 
[{^mt(k)}, {t/^m;(k)}, {t/'^ 4 ,(-k + Q)}, -{t/:;^^(-k + Q)}]*, 
with ma (k), dysCT (k), (k)]. The IRs 

for the onsite (k, —k + Q) pairings are summarized in 
Table III. Here the matrix for {*1^^}^ is 54 = S 0774 and 
^74 = 1 © —1 © (— 1 )“ with a = 1 for pairing and 

ct = — 1 for dy^-dxy pairing. The IRs for the non-onsite 
(k, —k + Q) pairings are summarized in Table IV. We can 
check that all the (k, —k + Q) pairing channels correspond 
to the inter-band pairings, and such kinds of pairings can not 
individually give a overall full gap around the Fermi surface. 


Appendix B: The classifications for the (k, —k + Q) pairing 
channels from symmetry analyses 


The eight GellMann matrices Aq Ag in the main text are 
listed as follows. 



■ 1 0 o' 


0 

1—1 

0 


Ao — 

0 1 0 

, ,^1 — 

1 0 0 

, ^2 — 


0 0 1 


0 0 0 




■ 1 0 o' 


'0 0 r 


'0 0 

^3 — 

0-10 

0 0 0 

, A 4 = 

0 0 0 

1 0 0 

, ,^5 — 

0 0 0 
i 0 0 


Ae — 

'0 0 o' 
0 0 1 

, ,^7 = 

'000' 

0 0 -f 


'100' 
0 1 0 


0 1 0 


0 i 0 

0 0-2 
(Bl) 
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